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1.  Introduction 

As  the  operating  frequency  of  the  millimeter- wave  circuits  is  increased,  some  of  the 
problems  that  are  not  very  relevant  at  lower  frequencies  become  significant.  For 
instance,  accurate  evaluations  of  the  discontinuities  in  the  planar  transmission  lines 
become  essential  in  the  design  of  the  integrated  circuits,  particularly  those  in  the 
monolithic  form.  It  is  no  longer  sufficient  to  simply  characterize  the  active  device  from 
the  DC  point  of  view.  Interactions  of  the  guided  waves  with  the  active  devices  becomes 
important  in  the  design  of  many  components.  It  is  desirable  and  necessary  to  explore 
new  types  of  circuits  and  structures  that  take  advantages  of  these  wave  interactions. 

The  projects  summarized  in  this  report  have  been  undertaken  with  the  spirit 
mentioned  above.  A  number  of  projects  on  the  interactions  of  the  guided  waves  with 
devices  have  been  studied.  This  report  summarizes  all  of  the  projects  supported  by  the 
Army  Research  Office  under  Contract  DAAG29-84-K-0076  and  the  publications  that 
resulted. 

2.  Planar  Transmission  Lines 

In  the  past  several  years,  numerous  publications  have  appeared  in  the  literature  on 
the  subject  of  guided  wave  characteristics  of  various  planar  transmission  lines.  In  spite 
of  all  this  activity,  accurate  characterizations  of  discontinuities  as  the  frequency  is 
increased  are  still  in  a  rather  primitive  stage  .  Additionally,  a  recent  discovery  has 
generated  some  interest  in  superconducting  planar  transmission  lines.  Projects  related  to 
these  subjects  will  be  summarized  below. 

2.1  Microstrip  Step  Discontinuities 

The  microstrip  step  discontinuity  problems  are  often  formulated  in  terms  of  the 
mode  matching  technique  applied  to  an  equivalent  waveguide  model.  Figure  1  shows  the 
step  discontinuity  as  analyzed  by  the  mode  matching  technique.  From  the  mode  matching 
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technique,  several  different  formulations  are  obtained  depending  on  how  the  E-  and  H- 
field  continuity  equations  are  processed.  Some  formulations  are  believed  to  be  superior 
to  other*  in  terms  of  convergence  and  numerical  labor.  An  extensive  study  has  been 
conducted  to  compare  all  possible  formulations  for  the  microstrip  discontinuity  problem. 
It  was  found  that  one  of  them  is  much  superior  to  all  others.  In  fact,  this  is  the  simplest 
formulation  and  requires  the  smallest  matrix  [1]  (see  Appendix  1).  To  study  the  accuracy 
of  the  method,  an  independent  check  of  the  results  has  been  conducted  by  using  the 
Modified  Residue  Calculus  Technique.  This  method  is  known  to  be  numerically  efficient 
and  to  result  in  accurate  scattering  parameters  of  the  discontinuity.  The  results  compare 
very  well  with  those  obtained  using  the  mode  matching  method  as  well  as  those  that  have 
been  reported  in  the  literature  [2]. 


Fig.l.  Microstrip  Step  Discontinuity 

When  the  mode  matching  method  or  the  Modified  Residue  Calculus  Technique  is 
combined  with  the  Generalized  Scattering  Matrix  Technique,  cascaded  and  offset  step 
discontinuities  can  be  characterized.  The  cascaded  step  discontinuities  often  appear  in  the 
form  of  a  symmetric  stub  to  the  main  microstrip  line  and  can  be  decomposed  to  two 
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successive  step  discontinuities.  On  the  other  hand,  the  offset  microstrip  discontinuity  can 
be  decomposed  to  two  successive  step  discontinuities  separated  by  a  zero  distance.  First, 
each  step  discontinuity  is  characterized  by  the  mode  matching  method  or  by  the  Modified 
Residue  Calculus  Technique  so  that  the  generalized  scattering  matrix  of  each  discontinuity 
is  derived.  These  generalized  scattering  matrices  are  combined  by  the  matrix 
manipulation  to  obtain  that  of  the  composite  structure.  The  interactions  between  the  two 
discontinuities  by  way  of  the  dominant  as  well  as  higher  order  modes  are  taken  into 
account  [3]  (see  Appendix  2). 


2.2  Non-Touching  E-Plane  Fin  and  Its  Application  to  an  Evanescent  Mode  Filter 
A  non-touching  E-plane  fin  shown  in  Fig.2  is  a  useful  circuit  component  for  E- 
plane  integrated  circuits.  It  is  printed  on  a  sustrate  which  is  in  turn  inserted  into  a  TEjq 

waveguide  along  its  E-plane.  One  end  of  the  stub  is  connected  to  the  top  or  bottom  wall 
of  the  waveguide  whereas  another  end  does  not  reach  the  opposite  wall.  This  element 
can  be  used  as  a  transition  from  a  waveguide  to  a  printed  circuit  and  as  a  filter  element. 
This  configuration  has  been  analyzed  by  an  extension  of  the  spectral  domain  technique. 
This  extension  essentially  provides  a  Fourier  transformed  version  of  the  integral 
equations  with  respect  to  unknown  current  components  excited  on  the  stub.  The 
formulations  are  algebraic  equations  with  inhomogeneous  terms  corresponding  to  the 
incident  mode.  These  equations  are  solved  for  the  current  components  on  the  stub.  Once 
these  current  components  are  found,  the  amplitude  of  a  particular  scattered  mode  can  be 
identified  from  the  residue  of  the  Green's  function  in  the  spectral  domain  [4,5]. 


Ill 
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Fig.2.  Non-Touching  E-Plane  Fin 


The  analysis  is  extended  to  the  case  where  a  non-touching  E-plane  fin  is  located  in 
the  waveguide  below  cutoff.  With  this  configuration,  an  evanescent  mode  filter  can  be 
formed.  The  evanescent  waveguide  containing  non-touching  E-plane  fins  are  connected 
to  the  standard  waveguides  on  both  ends  via  double  step  discontinuities  as  shown  in  Fig. 
3.  The  generalized  scattering  parameters  are  obtained  for  the  double  step  discontinuity  by 
means  of  the  mode  matching  method.  Each  fin  is  characterized  by  the  extended  spectral 
domain  method  and  the  generalized  scattering  parameters  are  obtained.  Connecting  all  of 
these  scattering  matrices,  one  can  find  the  scattering  matrix  of  the  entire  filter.  A  simple 
way  of  explaining  the  operation  of  this  filter  is  as  follows.  The  waveguide  below  cutoff 
is  essentially  a  quasi-lumped  inductance.  The  non-touching  E-plane  fin  provides  a 
capacitance  which  resonates  with  the  inductance  to  determine  the  passband  of  the  filter 
[6].  A  typical  passband  characteristic  is  shown  in  Fig.  4. 
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investigated  and  several  modifications  of  the  microstrip  line  have  been  suggested  for 
possible  dielectric  loss  reduction  [8]. 

3.  Quasi-Optical  Structures 

As  the  operating  frequencies  are  increased,  a  quasi-optical  structure  is  an  attractive 
alternative  for  conventional  microwave  and  millimeter-wave  systems.  A  typical  quasi- 
optical  system  consists  of  lenses,  reflectors  and  wire  grids.  The  wave  in  these  free- space 
waveguides  eventually  interface  with  microwave  circuits.  In  the  past  several  years,  there 
have  been  some  efforts  to  use  planar  circuit  technology  in  this  interfacing  portion  of  the 
structure.  A  number  of  quasi-optical  mixers  have  been  reported  in  which  a  planar 
antenna  is  directly  integrated  with  the  mixer  circuit  [9,10].  In  comparison  to  these 
mixers,  little  has  been  reported  on  studies  of  other  active  components.  Several  attempts 
have  been  made  in  this  area  during  this  contract  period. 


A  novel  structure  for  a  high-power  millimeter-wave  quasi-optical  application  has 
been  conceived  and  tested.  On  a  large  ground  plane,  a  number  of  periodically  placed 
slots  are  created  and  a  varactor  or  Schotkky  diode  is  placed  in  each  solt.  These  slots  are 
excited  by  a  reduced  height  waveguide  with  its  top  wall  as  a  part  of  the  ground  plane 
described  above  (see  Fig.  6.).  The  slots  are  resonant  at  the  desired  second  harmonic 
frequency  and  fed  by  the  fundamental  frequency  from  the  reduced  height  waveguide. 
Since  the  slots  are  electrically  small  at  the  fundamental  frequency,  only  a  small  amount  is 
coupled  to  each  slot.  Thus,  a  large  number  of  slots  can  be  fed  by  a  single  source.  No 
branching  network  or  tuning  circuits  are  required.  The  diode  can  be  operated  without  any 
dc  bias.  The  second  harmonic  signals  from  each  slot  are  radiated  into  free  space  and 
combined  there.  In  a  test  from  4  to  8  GHz  multiplication,  the  conversion  loss  is  about  2 
dB  and  in  the  one  from  17  to  34  GHz  it  is  7  dB.  Fourteen  diodes  are  used  in  the 
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structure  and  the  sidelobe  level  of  this  "phased  array"  is  less  than  -14  dB.  An  8-diode 
combiner  has  been  tested  for  doubling  a  35  GHz  signal  to  a  70  GHz  output  and  the 
sidelobe  level  was  less  than  -12  dB  [1 1]. 
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Fig.6.  Quasi-Optical  Multiplying  Slot  Combiner  Array 


3.2  Planar  Frequency-Doubling  Power  Combiner  Excited  by  a  Microstrip  Line 

The  frequency-doubling  power  combiner  is  fed  by  a  microstrip  in  this  project  so 
that  the  structure  is  more  amenable  to  planar  and,  eventually,  to  monolithic  fabrication. 

In  this  structure,  a  number  of  slots  in  the  ground  plane  on  the  substrate  are  excited  by  a 
microstrip  line  on  the  other  side  of  the  substrate  (see  Fig.  7).  First,  characteristics  of  the 
transition  from  a  microstrip  line  to  a  quarter  wave  slot  has  been  studied.  From  this  study, 
the  amount  of  power  left  over  to  the  output  side  of  the  microstrip  line  can  be  evaluated. 
The  location  of  the  transition  is  adjusted  at  each  slot  so  that  approximately  the  same 
amount  of  power  is  fed  to  each  slot.  Also,  the  length  of  the  microstrip  line  feed  between 
the  two  slots  is  adjusted  so  that  all  the  slots  are  fed  in  phase.  Based  on  these  studies,  slot 
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arrays  with  up  to  eight  elements  have  been  successively  tested.  The  conversion  loss  was 
about  16  dB  from  the  input  of  5.4  GHz  to  the  output  of  10.8  GHz  [12].  A  detailed 
design  consideration  is  reported  in  Appendix  3. 


(a) 


(b) 


Fig.7.  Microstrip  Fed  Quasi-Optical  Frequency  Doubling  Combiner 

3.3  Coupled  Slot  Antenna 

The  radiation  characteristics  of  a  pair  of  half  wave  slot  antennas  have  been  studied 
both  theoretically  and  experimentally.  This  project  is  to  confirm  that  the  even  mode 
(coupled  slot  mode)  radiates  well  while  the  odd  mode  (coplanar  wavegude  mode)  does 
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not.  The  calculations  are  based  on  the  spectral  domain  method  to  find  the  fields  in  the 
slot  as  well  as  in  the  far  field.  Experiments  of  the  radiation  patterns  agreed  well  with  the 
theoretical  prediction  [13]. 


Based  on  the  even  and  odd  mode  properties  of  a  coupled  slot,  a  self  oscillating, 
quasi-optical  mixer  has  been  designed  and  built  with  two  HEMT  devices  at  the  X  band. 
The  coupled  slot  receives  an  incoming  RF  signal  in  its  even  mode.  One  HEMT  is  placed 
in  each  slot  and  generates  a  local  oscillator  signal  which  is  supported  by  the  odd  mode  in 
a  coupled  slot.  The  RF  and  LO  signals  are  then  mixed  in  the  HEMT  devices  and  down 
converted  to  an  IF  signal  which  is  extracted  by  a  microstrip  line  on  the  substrate.  An 
isotropic  conversion  gain  of  5  dB  has  been  attained  [14].  Figure  8  shows  a  schematic  of 
the  structure. 
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Fig.8  Quasi-optical  HEMT  self  oscillating  mixer 
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4.  Device  Modeling  and  Characterizations 


At  higher  millimeter-wave  frequencies,  the  width  effect  of  the  three  terminal  devices 
such  as  MESFETs  can  be  significant  This  effect  may  be  detrimental  in  a  usual  device. 
However,  several  possibilities  exist  for  using  the  width  effect  to  create  a  distributed 
device  in  which  the  guided  wave  is  amplified  while  it  is  traveling  in  the  width  direction  of 
the  device.  Before  such  a  device  is  investigated,  it  is  necessary  to  characterize 
conventional  three  terminal  devices  both  from  the  device  simulation  point  of  view  and 
from  the  circuit  application  point  of  view. 

4.1  Large  Signal  Model  of  GaAs  MESFET 

An  efficient  numerical  algorithm  has  been  developed  for  a  large  signal  lumped 
circuit  model  of  the  GaAs  MESFET  which  will  aid  circuit  design  and  will  find  the  two 
port  charactersitics  of  the  device.  The  method  is  base  on  the  conventional  hybrid  time- 
frequency  domain  iterative  method  called  the  "multiple  reflection  technique."  However, 
the  convergence  has  been  improved  significantly  by  incorporating  an  acceleration 
algorithm  similar  to  the  successive  under-relaxation  technique.  The  predicted  results  for  a 
MESFET  power  amplifier  agree  well  with  the  measured  data  [15]. 

4.2  Two-Dipiensipnal  SimulatiaiLpf£ET  Structures 

FET's  of  several  different  types  have  been  simulated  by  means  of  a  finite  difference 
algorithm.  The  primary  purpose  of  this  work  is  to  select  a  candidate  for  a  traveling  wave 
FET.  The  simulation  program  takes  into  account  the  energy  momentum  relations  so  that 
the  method  can  simulate  FETs  with  a  submicron  gate  length  [16].  Among  these 
structures  investigated,  an  FET  with  an  inverted  gate  configuration  has  been  selected  for 
the  traveling  wave  device  (see  Fig.  9).  In  this  configuration,  the  gate  is  located  on  the 
underside  of  the  active  layer  through  the  substrate.  The  source  and  the  drain  are  located 
on  the  upper  side  of  the  active  layer  in  the  usual  manner.  The  device  should  be  operated 
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with  a  grounded  gate  mode.  Then  the  source  and  the  drain  can  be  treated  in  a  coupled 
strip  transmission  line  under  the  traveling  mode  of  operation.  Electrical  characteristics  of 
the  device  have  been  obtained  which  will  be  used  in  modeling  the  traveling  wave  device 
in  work  which  will  be  done  under  the  subsequent  contract. 
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5.  Slow  Wave  Structures 

One  of  the  problems  facing  the  implementation  of  millimeter-wave  monolithic 
integrated  circuits  is  the  large  size  of  passive  devices.  They  occupy  valuable  real  estate 
on  the  GaAs  substrate.  To  alleviate  this  problem  we  have  done  extensive  studies  on  the 
slow  wave  structures.  Two  types  of  slow  wave  structures  are  available.  One  of  them  is 
based  on  the  lossy  layer  in  a  semiconductor  substrate  and  another  makes  use  of  crosstie 
metal  configuration.  The  latter  one  is  proposed  for  extensive  study  in  the  subsequent 
contract.  Only  very  preliminary  work  was  done  on  this  during  the  report  period. 

5.1  Slow-Wave  Coupled  Line  Structures 

The  slow  wave  propagates  in  a  printed  transmission  line  in  which  the  conductors 
are  placed  on  a  thin  lossless  layer.  This  is  in  turn  placed  on  a  lossy  layer  which  is  created 
on  a  lossless  substrate  (see  Fig.  10).  The  cause  of  the  slow  wave  is  spatial  separation  of 
the  electric  energy  and  the  magnetic  energy.  The  electric  energy  is  mainly  concentrated  in 
the  thin  lossless  layer  while  the  magnetic  field  penetrates  into  the  lossy  layer  and  into  the 
substrate  beneath.  This  field  separation  can  be  made  by  an  appropriate  choice  of  the 
structural  and  material  parameters. 
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Fig.  10  Cross  Section  of  a  Slow  Wave  Structure 


The  slow  wave  propagation  in  a  coupled  line  can  be  controlled  by  a  number  of 
material  and  structural  parameters.  We  have  discovered  that  the  phase  velocities  of  the 
even  and  odd  modes  can  be  equalized  by  a  proper  choice  of  these  parameters.  Hence,  it 


is  possible  to  build  a  directional  coupler  with  very  small  physical  dimensions.  This 
feature  is  attractive  for  GaAs  monolithic  integrated  circuits  where  the  substrate  real  estate 
is  expensive.  A  computer  program  has  been  developed  for  design  of  directional  couplers 
both  in  the  forward  and  backward  waves  [17]. 

5.2  Crosstie  CPW  and  Microstrip  Slow  Wave  Structures 

A  new  slow  wave  structure  which  does  not  depend  on  the  substrate  loss  has  been 
conceived.  The  structure  is  either  a  coplanar  waveguide  (CPW)  or  a  microstrip  line  on 
which  a  thin  dielectric  overlay  is  placed.  Periodic  crosstie  strips  oriented  in  the  transverse 
direction  are  placed  on  the  overlay.  The  schematic  is  shown  in  Fig.  1 1.  The  capacitance 
between  the  crosstie  and  the  center  strip  of  the  microstrip  or  the  center  strip  of  the 
CPWcan  be  made  extremely  large.  Since  the  structure  is  now  a  cascaded  chain  of 
capacitive  and  inductive  sections,  the  electric  and  magnetic  energy  are  spatially  separated 
to  cause  a  slow  wave  effect.  It  is  important  that  the  operating  wavelength  be  sufficiently 
longer  than  one  period  of  the  crosstie  CPW  or  microstip  line  so  that  the  structure 
simulates  a  uniform  transmission  line.  Due  to  a  large  capacitance  and  inductance  per  unit 
length,  the  phase  velocity  is  extremely  small  and  hence  the  slow  wave  is  obtained.  A 
preliminary  calculation  suggests  that  a  slow  wave  factor  of  more  than  10  can  be  obtained 
[18].  There  are  a  number  of  conceivable  interesting  applications.  They  include  a 
bandstop  filter  made  of  a  doubly  periodic  grating  and  other  types  of  filters.  All  of  them 
are  expected  to  be  small  in  size.  These  devices  will  be  investigated  in  the  subsequent 
contract  period. 


Fig.  11.  Crosstie  Slow  Wave  Structrures 

6.  Broadband  Planar  Optical  Modulator 

The  major  limitation  of  a  traveling  wave  optical  modulator  is  velocity  matching 
between  the  guided  optical  wave  and  the  modulating  microwave  signal.  A  mismatch  of 
the  phase  velocity  between  these  two  waves  limit  the  bandwidth  of  the  modulator.  The 
optical  wave  is  guided  typically  by  a  waveguide  which  is  buried  in  a  substrate  material. 
On  the  other  hand,  the  microwave  signal  propagates  along  a  strip  conductor  on  the 
surface  of  the  substrate.  In  the  case  of  a  GaAs  based  modulator,  the  phase  velocities  of 
the  microwave  signal  is  much  faster.  Since  the  new  crosstie  overlay  coplanar  waveguide 
described  in  the  previous  section  has  the  capability  of  slowing  down  the  microwave 
signal,  this  structure  can  be  used  as  the  propagation  path  of  the  modulating  microwave 
signal  for  a  traveling  wave  integrated  optical  modulator  created  with  a  GaAs  material. 


Another  important  parameter  of  the  traveling  wave  modulator  is  the  charactersitics 
impedance  of  the  microwave  signal  path.  The  crosstie  structure  has  the  capability  of 
adjusting  the  characteristic  impedance  in  addition  to  the  phase  velocity.  It  is  possible  to 
attain  a  broadband  impedance  matching  by  this  structure. 

The  crosstie  structure  is  provided  via  a  dielectric  overlay  on  a  coplanar  waveguide 
electrode  which  is  placed  on  a  GaAs/GaAlAs  layered  structure.  The  optical  waveguide  is 
in  the  form  of  a  GaAs  channel  waveguide  (see  Fig.  12).  A  preliminary  investigation 
leading  to  the  design  of  this  optical  modulator  is  reported  in  [19]  (see  Appendix  4). 
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Fig.  12.  Crosstie  CPW  Traveling  Wave  Optical  Modulator 
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Comparative  Study  of  Mode-Matching 
Formulations  for  Microstrip 
Discontinuity  Problems 
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Abstract  —  Several  matrix  formulations  for  the  microstrip  step-discon¬ 
tinuity  problem  are  compared.  Although  they  are  theoretically  identical, 
one  of  them  has  an  advantage  in  numerical  labor,  relative,  and  absolute 
convergence.  Results  of  this  method  are  checked  with  other  published  data 
and  with  those  independently  obtained  by  the  modified  residue  calculus 
technique. 

I.  Introduction 

A  STEP  DISCONTINUITY  is  frequently  encountered 
in  microstrip  line  circuits  and,  hence,  its  analysis  is 
important  for  circuit  design.  There  are  several  approaches 
available.  When  the  microstrip  is  enclosed  in  a  waveguide¬ 
like  case,  it  is  possible  to  calculate  the  fundamental  and 
higher  order  modes  in  both  sides  of  the  step  and  to  impose 
the  continuity  conditions  of  the  tangential  field  in  the  cross 
section  of  the  shield  case  at  the  step  location.  This  process 
leads  to  a  system  of  mode-matching  equations.  When  an 
open  microstrip  line  circuit  is  dealt  with,  the  higher  order 
modes  can  become  radiation  modes  which  must  be  in¬ 
cluded  in  the  mode-matching  procedure. 

In  many  applications,  the  so-called  waveguide  model  has 
been  found  useful  for  calculation  of  the  scattering  at  the 
microstrip  discontinuity  (1J— [3].  In  this  paper,  we  assume 
that  the  waveguide  model  is  an  acceptable  technique  for 
the  calculation  of  the  step-continuity  problem.  In  addition, 
radiation  and  surface  waves  are  not  considered.  The  moti¬ 
vation  for  the  present  work  is  somewhat  different  from 
those  published.  Although  a  number  of  numerical  data  are 
presented  in  the  literature  [1  ]— [3J,  details  of  the  numerical 
process  are  not  clear.  The  objective  of  the  present  paper  is 
not  to  duplicate  the  numerical  data  already  available,  but 
to  place  some  foundation  on  how  these  data  should  be 
calculated.  We  present  several  alternative  formulations. 
Although  these  formulations  are  theoretically  identical,  it  is 
pointed  out  that  the  numerical  labor  and  accuracy  depend 
on  the  choice  of  formulation  and  some  are  better  than 
others. 
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Fig  1.  The  waveguide  model  of  the  open  microstrip  line. 

The  best  formulation  can  be  decided  based  on  the  matrix 
size,  relative  and  absolute  convergence,  and  other  numeri¬ 
cal  considerations.  It  turns  out  that  the  best  formulation  is 
the  one  we  often  choose  without  clear  reasoning.  The  data 
for  a  microstrip  step  discontinuity  are  compared  with 
available  data.  They  are  also  compared  with  the  modified 
residue  calculus  technique,  which  serves  as  an  independent 
check  of  the  numerical  accuracy. 

Before  starting  the  formulation,  let  us  briefly  review  the 
waveguide  model.  In  this  technique,  the  uniform  microstrip 
line  of  width  w0  on  the  substrate  of  height  h  and  relative 
dielectric  constant  tr  is  replaced  with  an  equivalent  paral¬ 
lel-plate  waveguide  with  magnetic  side  walls  (Fig.  1).  The 
substrate  height  is  kept  identical.  However,  the  effective 
width  weff  and  the  effective  dielectric  constant  teff  are  used 
to  define  the  effective  waveguide  in  such  a  way  that  the 
effect  of  the  fringing  field  of  the  microstrip  is  taken  into 
account.  Specifically,  these  effective  values  are  related  to 
the  propagation  constant  fi  and  the  characteristic  imped¬ 
ance  Z0  via 

€eff  “  ( P/ko) 

Z0  =  (l20ir/Q(  */*,„)•  (1) 

Note  that  /3  and  Z0  can  be  found  from  a  standard  analysis 
such  as  the  spectral-domain  method,  from  curve  fitting,  or 
an  empirical  formula,  once  the  structural  parameters  of  the 
microstrip  line  are  given. 

For  the  analysis  of  the  step  discontinuity,  both  sides  of 
the  step  are  replaced  with  their  respective  equivalent  wave¬ 
guides.  Note  that  the  heights  of  these  two  waveguides  are 
identical.  Hence,  the  problem  remains  a  two-dimensional 
one  as  the  field  is  uniform  in  the  y-  (vertical)  direction. 
Also,  note  that  the  dominant  mode  in  the  equivalent 
waveguide  is  TEM.  In  the  following  sections,  all  structural 
parameters  used,  except  h,  are  presumed  to  be  the  “effec¬ 
tive”  ones,  unless  otherwise  stated. 
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Fig.  2.  (a)  Waveguide  model  for  symmetric  microstrip  step  discontinu¬ 
ity.  (b)  Auxiliary  geometry  for  the  waveguide  problem. 


II.  Formulation 

The  problem  under  study  is  the  waveguide  model  for  the 
microstrip  step  discontinuity  shown  in  Fig.  2(a).  The  struc¬ 
ture  is  assumed  to  be  symmetrical,  and  the  parallel-plate 
waveguide  is  idealized  with  magnetic  side  walls.  For  con¬ 
venience  of  analysis,  an  auxiliary  structure  is  introduced  in 
Fig.  2(b).  Only  one  half  of  the  original  structure  is  consid¬ 
ered  because  of  symmetry,  and  the  transversal  magnetic 
wall  at  the  discontinuity  is  recessed  to  create  a  new  region 
C.  The  original  structure  is  recovered  by  letting  d  =  0. 

The  mode-matching  procedure  begins  by  expanding  the 
tangential  electric  and  magnetic  fields  at  the  junction  in 
terms  of  the  normal  modes  on  both  sides  of  the  junction. 
For  TE„0  (n  =  0,1,  •  •  ■ )  excitation,  we  write  down  the  Ey 
continuity  equation 


L  (An  +  A~  Ha*  ”  '  i-1 

""°  £CA„(1  +  P«).  bZx<,a 

,  it  -  o 


and  a  corresponding  one  for  Hx 


where 

4>o*-A(*o/a)c°s(kamx),  kan  **  (rtir/a) 

<bbn~‘A(^o/b)cos(kb„x),  kbm-{nv/b) 

<f>c-v(<,.oA)co  >s(kCH(a-x)),  kCH  =  (rt-rr/c) 

«„0-l,  n-0 

2,  n  +  0 

and 

Ybn  =  Af2kl-(™/b)1)  =  yn 

rcnxA(iko-(nir/c)2)  =  /?„ 

ft,~«p(-2y/U). 

In  (2),  4>bn,  and  <f>cn  are  normal  modes  in  Regions  A, 

B,  and  C,  with  propagation  constants  /?„,  y„,  and  /}„, 
respectively.  A*  and  B~  are  the  given  incident  field  coeffi¬ 
cients  from  Regions  A  and  B  (usually  only  one  A*  or  B~  is 
considered  at  a  time),  while  A~,  B*,  and  C„  are  the 
unknown  excited  field  coefficients  in  regions  A,  B,  and  C, 
respectively.  pn  is  the  reflection  from  the  magnetic  wall  in 
Region  C. 

From  modal  orthogonality,  we  obtain  the  linear  simulta¬ 
neous  equations  for  the  unknown  modal  coefficients 

a: + a- = K +b„-)+ l£hm\ + pj 

*  —  0  n  -  0 

Yam(A : -a~)  =  K£ hMb;  -  B; ) 

n  —  0 

+  E 

m  =  0,1,2, •  •  • ,  M  —  1  (3a) 

L  Hnm{A:  +  A;)  =  B:  +  Bm 

n  —  0 

t 1  Hmmr.*(A;  -a;)-  Yhm(  b:  -B-), 


m  =  0,1,2,-  •  • ,  A  - 1  (3b) 


L  #,m(J:  +  A;)-cm(i  +  pJ 


L  H„mYajA;-A;)-CmYcm(\-pm), 


m  =  0,1,2,  ■  •  • ,  L  - 1  (3c) 


where 


Z  (B;-B„-)Ybn<pbn,  o  sxsi 


E  CU.(l‘fl.),  bzxza  (2b) 


!" ^ /V •  >jv>> ;/  vV; >;>v ^ 


j\m+h*dx 

■'o 

_  j1/2|  1  2/Cam(~l)n$in(/c^Z))  \ 

\  ^  aw  ^  I 

_  /-a 

H  —  <t)  d>  dx 

aImn  I 
Jh 

_  /  <motnO\1/2(  1  lkam(-\)'*\m(kanb 
>  '2'  k2am  -  k2cn 
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To  condense  the  above  equations,  we  define  the  following 
matrices: 


0 

y,x 


1  +  p0  0 
0  l  +  pt 


i-a,b,c 


1-Po 


/?'  = 


0 


0 

1-Pi 


1  +  P„ 


1-PBJ 

n  »  L  —  1 


Y  1 
*b  I 

0 

i ;  o 

7 ;  o 

1 

- 1 

! 

[h  ! 

~n. 

H } 

R‘  = 

R  = 

oT  R 

where  /  is  the  identity  matrix,  H  is  a  matrix  of  size  3/  x  K 
with  generic  element  Hm„  as  defined  above^while  H  is  a 
matrix  of  size  M  x  L  with  generic  element  Hmn. 

Then,  the  mode-matching  equations  can  be  written  in 
the  following  matrix  form: 

a*  +  a~  =  GRd*  +  Gd~  (4a) 


Ya{q*  -q-)~GYdR'd+  -  GYdd~  (4b) 


Gr(a*  +  a~  )  *•  Rd*  +  d~  (4c) 

GTYa(q*-q-)-YdR’d/-Ydd-  (4d) 


M 


25 

Fig.  3.  Oassification  of  formulations. 


B0+ 


a*  and  d~  are  column  vectors  of  the  excitation  terms  and 
a~  and  d *  are  column  vectors  of  the  unknown  modal 
coefficients.  All  matrices  are  of  size  (A/  x  A/);  this  requires 
that  K  +  L~M. 

When  M  -»  oo,  we  can  prove  that  (7~ 1  =  Gr.  Therefore, 
(4a)  and  (4b)  are  equivalent  to  (4c)  and  (4d).  Two  indepen¬ 
dent  vector  equations  are  required  to  solve  for  two  un¬ 
known  vectors.  Hence,  for  four  pairs  of  equations  ((4a)  and 
(4b),  (4b)  and  (4c),  (4c)  and  (4d),  and  (4d)  and  (4a)). 
substituting  one  equation  into  the  other  in  the  same  pair, 
we  have  eight  ways  to  solve  for  a*  and  d~.  They  are 
defined  graphically  in  Fig.  3.  The  approaches  indicated  by 
a  solid  arrow  are  classified  as  the  formulations  of  the  first 
kind  and  those  indicated  by  a  dashed  arrow  are  of  the 
second  kind.  Although  the  eight  ways  of  solution  are 
theoretically  equivalent,  their  numerical  behaviors  are 
somewhat  different,  especially  when  the  magnetic  wall  is 
introduced  at  the  upper  half  of  the  junction  (d  =  0.  p„  =  1). 

For  general  cases  d  *  0,  all  the  formulations  require  a 
matrix  inversion  of  size  (M  X  M).  For  our  limiting  case  of 
d  =  0,  special  modifications  must  be  made  for  some  cases. 
Specifically,  ID  and  2B  need  to  invert  a(M  +  L)x(M  +  L) 
matrix  and  2C  (Appendix  A)  needs  to  invert  a  smaller 
(K  X  K)  matrix.  Hence,  2C  is  most  attractive  to  us  be¬ 
cause  of  its  potential  of  numerical  efficiency.  In  the  next 
section,  we  will  examine  the  various  approaches  in  terms  of 
the  numerical  stability  and  convergence. 

III.  Numerical  Results 


where  superscript  T  denotes  transpose  operation,  and 


< 

’  *o  ' 

a+  — 

d '  = 

Bk- i 

"  "o" 

To  study  the  numerical  behavior  of  the  various  formula¬ 
tions,  we  have  chosen  the  structural  parameters  as:  a  = 
100,  b  =  26.1  (in  mils),  <,  =  2.2, t2  =  2.1.  The  dominant 
mode  (TEM)  reflection  and  transmission  coefficients  at  the 
junction  are  calculated  by  varying  the  matrix  size  for 
different  K/M  ratios. 

Since  formulations  ID  and  2B  have  an  apparent  disad¬ 
vantage  in  numerical  calculations,  they  are  not  considered 
here.  After  extensive  studies,  we  have  found  that  1A,  IB. 
and  1C  are  numerically  identical.  Similarly,  2A  and  2D  are 
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Fig.  4.  Convergence  study  for  various  formulations:  (a)  formulation  1A; 
(b)  formulation  2A;  (c)  formulation  2C 

numerically  identical.  Therefore,  only  three  sets  of  data, 
corresponding  to  1A,  2A,  and  2C,  are  given. 

In  each  formulation,  the  indices  L ,  K,  and  M  are 
involved.  The  numerical  results  are  affected  by  the  ratios 
among  these  indices.  This  is  called  the  relative  convergence 
phenomenon,  and  it  has  been  thoroughly  discussed  in  the 
literature  (4],  [5].  It  is  well  known  that  the  best  approxima¬ 
tion  to  the  true  solution  is  obtained  for  L/K  =  c/b  = 
73  6/26.1  or  M/K  ■=  a/b  =  100/26.1  (refer  to  as  the  right 
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(a) 


(b) 

Fig.  5  Relative  convergence  problem  of  formulation  1A  demonstrated 
by  field  plots. 


ratio).  It  is  observed  that  2A  and  2C  suffer  very  little  from 
the  relative  convergence  problem  (Fig.  4(b)  and  (c)).  The 
problem  is  more  serious  in  1A,  as  can  be  seen  in  Fig.  4(a). 
With  a  ratio  of  one,  the  dominant  mode  reflection  coeffi¬ 
cient  converges  to  an  incorrect  value  (curve  A  in  Fig.  4(a)). 
Curve  C  is  calculated  using  a  ratio  of  three,  which  is  close 
to  the  right  ratio. 

The  relative  convergence  effect  can  be  more  clearly 
observed  from  the  plot  of  Hx  at  the  junction.  The  resultant 
field  calculated  by  1A  is  shown  in  Fig.  5.  In  Fig.  5(a),  we 
plot  the  fields  calculated  using  a  ratio  of  L/K  =  28/10 
(very  close  to  the  right  ratio)  and  a  ratio  of  29/9  (higher 
than  the  right  ratio).  In  Fig.  5(b),  we  compare  the  field 
calculated  using  a  ratio  of  27/11  (lower  than  the  right 
ratio)  and  a  ratio  of  28/10.  It  is  interesting  to  note  that, 
with  a  higher  ratio,  the  calculated  field  behaves  better  on 
the  magnetic  wall  discontinuity  than  on  the  aperture,  while 
the  opposite  is  true  for  field  calculated  using  a  lower  ratio. 
This  might  seem  reasonable  because,  for  a  higher  ratio,  we 
use  more  modes  on  the  magnetic  wall  than  on  the  aperture 
and  vice  versa  for  a  lower  ratio.  Fig.  6  shows  the  resultant 
fields  calculated  by  2C.  Different  ratios  have  no  noticeable 
effect.  The  fields  calculated  by  2A  behave  similar  to  those 
calculated  by  2C. 
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Fig.  6.  Plot  of  fields  calculated  using  formulation  2C  with  different 
M/K  ratios. 


Fig.  7.  Comparison  of  numerical  efficiency. 

A  comparative  study  on  the  numerical  efficiency  for 
different  approaches  has  also  been  done.  In  this  case, 
M/K  =  4,  which  is  close  to  a/b,  is  chosen.  The  results  of 
the  dominant  mode  reflection  coefficient  S[00j(00]  are 
evaluated  as  a  function  of  the  matrix  size  required  and 
shown  in  Fig.  7.  In  addition,  a  comparison  of  how  well  the 
fields  of  the  two  sides  match  at  the  junction  is  done 
between  formulations  2A  and  2C.  In  both  calculations,  M 
is  set  to  40,  K  to  10.  The  result  is  shown  in  Fig.  8.  The 
fields  calculated  by  2C  match  as  well  as,  if  not  better  than, 
those  calculated  by  2A.  Keep  in  mind  that  we  have  to 
invert  a  matrix  of  size  40  in  2A,  compared  to  a  matrix  of 
size  10  in  2C.  It  is  now  obvious  that  2C  has  definite 
advantages  over  other  approaches.  This  formulation  is  to 
be  used  for  further  studies. 

Let  us  refer  back  to  (2)  at  this  point.  In  many  attempts, 
Ey  in  the  region  b  <  x  <  a  is  not  used  as  Hx  =  0  there.  This 
is  identical  to  2C  and,  hence,  is  our  preferred  choice. 

To  check  the  validity  of  our  calculations,  we  have  calcu¬ 
lated  the  frequency  response  of  a  microstrip  step  discon- 


(b) 

Fig.  8.  Comparison  of  field  plot  between  formulation  2A  and  2C.  (a) 
2A.  (b)  2C. 


tinuity  using  the  same  parameters  given  by  Kompa  [2],  The 
results  of  the  dominant  mode  reflection  and  transmission 
coefficients  S(00][00]  and  S[00][00]  are  shown  in  Fig.  9 
They  are  in  good  agreement  with  Kompa’s  results.  The 
small  discrepancy  is  due  to  the  different  formulas  used  for 
obtaining  the  effective  width  and  effective  dielectric  con¬ 
stant  of  the  waveguide  model.  Furthermore,  we  have 
checked  the  results  with  those  independently  obtained  by 
the  modified  residue  calculus  technique  [6].  The  results  of 
S[00][00]  are  shown  in  Table  I  and  Fig.  10  for  comparison. 
The  calculations  are  performed  using  Kcmpa's  parameters. 


IV.  Conclusion 

The  mode-matching  method  has  been  applied  to  analyze 
the  microstrip  step-discontinuity  problem  based  on  the 
waveguide  model.  A  comparison  has  been  made  among  the 
various  mode-matching  solutions  based  on  the  matrix  size 
and  relative  and  absolute  convergence.  Although  they  are 
theoretically  identical,  one  of  them  proves  to  be  most 
suitable  for  numerical  calculations.  The  results  by  this 
method  are  in  good  agreement  with  other  published  data 
and  with  those  independently  obtained  by  the  modified 
residue  calculus  technique. 
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Fig.  9.  Comparison  with  Kompa’s  results. 
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Fig.  10  Comparison  with  results  of  modified  residue  calculus  tech 
tuque. 


TABLE  f 

Comparison  of  the  Results  by  Mode-Matching  Method  and 
by  MRCT* 


VoOKlio) 


Mode  Matching 
01837-  /0.02291 
-0.7856+  j  0.2121 


MRCT 

0.1837-  ;0. 02297 
-0.7855  +  y0.2233 


R.  Mittra,  T.  Itoh,  and  T.  S.  Li,  “Analytical  and  numerical  studies  of 
the  relative  convergence  phenomenon  arising  in  the  solution  of  an 
integral  equation  by  the  moment  method,”  IEEE  Trans.  Microwave 
Theory  Tech.,  vol  MTT-20,  pp.  96-104,  Feb.  1972. 

T.  S.  Chu  and  T.  Itoh,  “Analysis  of  microstrip  step  discontinuity  by 
the  modified  residue  calculus  technique,”  submitted  to  IEEE  Trans. 
Microwave  Theory  Tech.,  (Special  Issue  on  Numerical  Technique). 
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Generalized  Scattering  Matrix  Method  for 
Analysis  of  Cascaded  and  Offset  Microstrip 

Step  Discontinuities 

TAK  SUM  CHU  and  TATSUO  ITOH,  fellow,  ieee 


Abstract  —  Detailed  algorithms  are  presented  for  characterizations  of 
cascaded  microstrip  step  discontinuities,  symmetric  stubs,  and  offset  step. 
The  analysis  is  based  on  the  generalized  scattering  matrix  techniques  after 
the  equivalent  waveguide  model  is  introduced  for  the  microstrip  line. 

I.  Introduction 

HE  EQUIVALENT  waveguide  model  has  been  ad¬ 
vantageously  used  for  characterizing  a  number  of 
discontinuities  appearing  in  microstrip  circuits  (l]-[3). 
Although  the  radiation  and  surface-wave  excitation  are 
neglected,  these  characterizations  provide  useful  and  accu¬ 
rate  information  in  many  practical  applications.  Formula¬ 
tions  for  these  discontinuity  problems  are  typically  done  by 
the  mode-matching  technique.  However,  detailed  formula¬ 
tion  algorithms  are  not  readily  available. 

Recently,  the  present  authors  carried  out  an  assessment 
for  a  number  of  different  formulations  for  a  microstrip 
step  discontinuity  within  the  framework  of  the  waveguide 
model.  The  most  economical  and  yet  most  accurate  formu¬ 
lation  was  suggested  [4].  The  step  discontinuity  has  also 
been  analyzed  by  the  modified  residue  calculus  technique 
(MRCT)  [5], 

The  present  paper  extends  the  analysis  of  the  step  dis¬ 
continuity  to  a  cascaded  step  discontinuity  and  an  offset 
step  discontinuity.  The  symmetric  stub  can  be  treated  as 
the  cascaded  step  discontinuity.  The  offset  discontinuity 
will  be  treated  as  the  limiting  case  of  a  cascaded  discon¬ 
tinuity.  In  each  case,  the  individual  step  discontinuity  is 
characterized  by  either  the  MRCT  or  the  mode-matching 
method.  The  analysis  results  in  a  generalized  scattering 
matrix  for  each  step.  The  analysis  of  the  cascaded  step  is 
undertaken  by  invoking  the  generalized  scattering  matrix 
technique  in  which  the  generalized  scattering  matrices  of 
two  step  junctions  are  combined  [6],  [7],  The  ultimate  result 
is  the  generalized  scattering  matrix  of  the  cascaded  junc¬ 
tion  as  a  whole. 

It  should  be  noted  that  the  waveguide  model  is  presumed 
to  be  an  acceptable  model  for  the  present  analysis.  The 
results  are  compared  with  the  experimental  data  as  well  as 
those  reported  in  the  literature. 

Manuscript  received  August  2,  1985;  revised  October  1,  1985  This 
work  was  supported  in  part  by  the  U  S  Army  Research  Office,  under 
Contract  DAAG29-84-K-0076 
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II.  Formulation  of  the  Problem 

A.  Cascaded  S.ep  Discontinuity 

The  algorithm  for  the  analysis  will  be  best  iilustr 
means  of  the  cascaded  step  discontinuity.  Specific  i 
required  for  the  offset  discontinuity  will  be  explainc 
The  first  step  of  the  analysis  is  to  replace  the  mil 
circuit  under  study  with  its  equivalent  waveguide  m< 
The  top  and  bottom  are  electric  walls  and  the  sidew 
magnetic  walls.  The  height  h  in  Fig.  1  remains  unci 
The  effective  dielectric  constant  «!  and  the  effectiv 
Wx  of  region  I  with  the  microstrip  width  can  b< 
from  the  dominant  mode  phase  constant  /?,  and  th 
acteristic  impedance  Z01  of  the  microsinp  line 

*i  =  (P  i/^o) 

Z01=[l20*/y'77](V^). 

and  Zoi  must  be  calculated  from  the  structural  p 
ters  by  a  standard  full-wave  analysis  [9],  [10)  or  a  c 
formula  [11],  Other  regions  may  be  modeled  in  a 
manner. 

The  next  step  is  to  characterize  all  the  discont 
involved  in  the  waveguide  model  of  the  microstnp 
under  study.  This  characterization  is  done  in  term; 
generalized  scattering  matrix  [6],  [7],  This  matrix  is 
related  to  the  scattering  matrix  used  in  the  mic 
network  theory,  but  differes  in  that  the  dominant  as 
higher  order  modes  are  included  Therefore,  the 
alized  scattering  matrix  will  be.  in  general,  of  infiniti 
Consider,  for  instance,  the  TE^q  excitation  with  unit 
tude  from  the  left  to  junction  1  in  Fig.  2.  If  the  am 
of  the  nth  mode  of  the  reflected  wave  to  the  left  is 
( n,p )  entry  of  the  scattering  matrix  Sn(n,p) 
Similarly,  if  the  amplitude  of  the  mth  mode  of  th 
transmitted  to  the  right  in  B„,  S2l(m,  p)  is  Bm 
matrix  elements  can  be  derived  similarly  Hence,  th< 
alized  scattering  matrix  S,  of  the  junction  1  can  be 
in  terms  of  four  submatrices  of  infinite  order 


The  corresponding  matrix  of  the  junction  2  is 

c  fs”  SU1 
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Fig,  l  Cascaded  step  discontinuities,  equivalent  waveguide  model,  and 

top  view 


All  of  the  elements  of  the  generalized  scattering  matrix  can 
be  obtained  by  solving  the  electromagnetic  problem  of  the 
junction  scattering  by  means  of  a  standard  technique  such 
as  the  mode-matching  technique  [4]  or  a  modified  residue 
calculus  technique  [5],  Since  the  details  of  these  methods 
for  an  isolated  step  discontinuity  are  reported  in  [4]  and 
[5J,  they  are  not  repeated  here.  We  presume  all  of  these 
quantities  are  now  available. 

The  remaining  step  is  to  combine  these  generalized 
scattering  matrices  of  the  cascaded  junctions  and  to  derive 
the  composite  matrix. 


SAA  SAC 
SCA  Scc 


(5) 


To  this  end,  we  introduce  the  transmission  matrix  SU  )  of 
the  waveguide  between  junctions  1  and  2.  The  wave  travels 
a  distance  L  so  that  each  mode  is  multiplied  by  exp{  -  y L) 


£</.)  = 


.-Vli 


Y.  1 


(6) 


where  y„  is  the  propagation  constant  of  the  nth  mode  of 
Region  B.  Hence  y,  =  jf}2  L ,  where  fi2  is  the  dominant-mode 
phase  constant  of  Region  B.  Our  algebraic  process  to 
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Fig.  2.  Derivation  of  S'  parameters  for  the  cascaded  step  discontinuity 
structures 


derive  A  is  detailed  in  the  Appendix.  The  results  are 


SAA  ~  s11  +  SnSUAU2S"S{l  'S1' 

(7a) 

SAC  =  S12S(,,l/2S34 

(7b) 

SC4  =  S43S,'',t/]S21 

(7c) 

SCC  =  544  +  S41S,  nU[S:2Si  1  >5  ’4 

(7d) 

where 

t/1  =  (/-S22S<uS33S,nrl 

(8a) 

U2  =  (1-S”SU)SUS{I))  1 

(8b) 

and  /  is  the  unit  matrix.  The  above  matrices  are  formally 
of  infinite  size.  However,  in  practice,  these  matrices  must 
be  truncated  to  a  finite  size.  It  is  found  that  excellent 
convergence  is  obtained  when  3x3  or  even  2x2  sub- 
matrices  S".  etc.,  are  used. 

It  should  be  noted  here  that  the  use  of  generalized 

scattering  matrices  is  increasingly  more  important  as  the 
distance  between  two  junctions  is  smaller  Therefore,  the 
present  technique  can  be  used  for  analysis  of  the  symmetric 
stub  from  the  knowledge  of  the  generalized  scattering 
matrix  of  a  single  step  discontinuity. 

B  Offset  Step  Discontinuity 

Next,  the  technique  described  above  will  be  applied  to 
an  offset  step  discontinuity  shown  in  Fig.  3.  This  offset 
discontinuity  occurs  in  a  microstrip  ctrcuit  either  intention¬ 
ally  or  unintentionally.  As  we  will  see  shortly,  a  small 
amount  of  offset  A  significantly  affects  the  scattering 
characteristic  of  the  discontinuity. 
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(b) 


Fig.  4  Equivalent  waveguide  model  of  the  offset  step,  (a)  perspective 
view  and  (b)  top  view 


Once  again,  the  first  step  is  to  derive  the  equivalent 
waveguide  model  depicted  in  Fig.  4.  Although  a  direct 
formulation  of  Fig.  4  is  possible  by  way  of  the  mode¬ 
matching  technique,  we  will  take  an  alternative  approach 
in  which  the  formulations  and  the  generalized  scattering 
matrices  of  the  symmetric  step  discontinuities  are  advanta¬ 
geously  used.  To  this  end,  an  auxiliary  structure  in  Fig.  5  is 
introduced.  Notice  that  the  original  offset  step  discontinu¬ 
ity  structure  can  be  recovered  by  letting  5  in  Fig.  5  to  zero 
after  all  the  formulations  are  carried  out.  Also,  the  individ¬ 
ual  discontinuities  71  and  72  in  Fig.  5  are  one-half  of  the 
symmetric  step  discontinuities.  Hence,  all  of  the  previous 
results  for  the  symmetric  step  discontinuities  excited  by  the 
even-mode  can  be  directly  used.  In  fact,  in  [4]  and  [5],  only 
one-half  of  the  structure  has  been  used  for  analysis. 

Once  the  scattering  matrices  of  71  and  72  are  available, 
the  scattering  matrix  of  the  composite  discontinuity  can  be 
derived  from  (7)  and  (8)  except  that  S{L)  =  /  when  6  — •  0. 
This  completes  the  formulation  for  the  offset  discontinuity. 
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Fig.  6.  Numerical  data  for  the  cascaded  step. 


III.  Results  and  Discussions 

,  Fig.  6  shows  typical  results  of  the  amplitude  of  the 
dominant  mode  (quasi-TEM)  reflection  coefficient  from  a 
cascaded  step  discontinuity.  The  results  are  found  to  agree 
very  well  with  the  experimental  data  taken  at  Hughes 
Torrance  Research  Center  for  the  microstrip  circuit  on  a 
Duroid  substrate. 

For  efficient  calculation,  it  is  desirable  to  be  able  to 
truncate  the  matrices  at  as  small  a  size  as  possible  with 
accurate  results.  In  the  case  of  the  cascaded  discontinuity, 
SiL)  contains  a  convergent  factor  since  all  of  the  higher 
order  modes  have  real  values  of  y„.  In  the  case  of  an  offset 
discontinuity,  such  exponentially  decaying  factors  disap¬ 
pear  because  the  length  8  -*  0,  and  hence  SU  )  =  /.  To  test 
the  convergence  of  the  solution,  the  dominant-mode  trans¬ 
mission  coefficient  calculated  using  generalized  scattering 
matrices  of  sizes  2x2,  3x3,  and  5x5  are  compared. 
Physical  parameters  are  chosen  to  be  the  same  as  those 
studied  by  Kompa  [12]  so  as  to  permit  a  comparison  of 
results.  Fig.  7  shows  the  results  of  this  convergence  study. 
It  is  seen  that  even  the  2x2  matrix  gives  reasonably 
accurate  results.  To  establish  the  validity  of  the  results, 
they  are  compared  with  those  calculated  by  Kompa  [12]. 
This  comparison  is  shown  in  Fig.  8(a)  and  (b).  Finally,  the 
dominant-mode  transmission  coefficient  for  various  ec¬ 
centricities  are  calculated.  It  is  evident  in  Fig.  9  that  the 
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Fig  7  Convergence  study  of  the  generalized  scattering  matrix  technique 
for  an  offset  step:  (a)  magnitude  and  (b>  phase. 
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Fie.  8  Comparison  with  Kompa's  results:  (a)  magnitude  and  (b)  phase. 


Fig.  9  Effect  of  eccentricity:  (a)  magnitude  and  (b)  phast 


effect  of  the  first  odd-mode  cutoff  is  exhibited  as  soon  as 
the  eccentricity  is  nonzero.  Also,  for  A  =  0.5  cm,  the  effect 
of  the  second-  (even-)  mode  cutoff  is  quite  small  due  to  the 
fact  that  the  center  of  the  smaller  waveguide  is  located  at 
the  second-mode  null  and,  hence,  little  coupling  via  this 
mode  exists. 


IV.  Conclusions 

The  generalized  scattering  matrix  technique  has  been 
applied  to  the  problems  of  a  cascaded  step  discontinuity 
and  an  offset  step  discontinuity.  The  waveguide  model  has 
been  presumed  to  be  applicable  for  analysis.  Individual 
discontinuities  are  characterized  first  and  the  results  are 
used  for  the  description  of  the  composite  discontinuity  via 
the  generalized  scattering  matrix  technique.  In  spite  of  the 
two-step  process,  the  overall  effort  in  numerical  processing 
is  quite  efficient.  In  actual  numerical  software,  the  pro¬ 
grams  for  isolated  discontinuities  can  be  used  as  sub¬ 
routines  for  the  composite  structures. 
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Appendix 

Derivation  of  Scattering  Matrix  for  the 
Cascaded  Step  Discontinuities 

The  following  matrix  equations  can  be  derived  from  Fig. 


’  RiA) 

\s11 

s12' 

J{A) 

Ra) 

is21 

s22. 

.  I'2) . 

0  Sa,]\ 

<n  o  JL 


R0)  _ 

S33 

s34 

'  ;<3) 

R(C) 

s43 

s**. 

Next.  (A2)  is  substituted  into  (Al)  to  get 

R(A)  =  S,ll{/4)  +  Sl2Sil)R°l  (A4) 

RM  =  S2lI'A)  +  S22S{l)R0).  (A5) 

Next,  (A2)  is  substituted  into  (A3)  to  get 

R0)  =  S33S(,7?(2'  +  S34/(°  (A6) 

R<C)  =  S4}SiL)R{2)  +  S44IU).  (A7) 

Equations  (A5)  and  (A6)  are  used  to  isolate  Rl2)  and  R0) 
Ra)  =  l^S21/'4'  +  l/,S22S(,)S34/<n  (A8) 


where 


Ul^(I~S22S,l)S}yS,u) 


R°'  =  U,SnS{l  IS21/*4’  +  U,SM1'C) 


where 
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i/:  =  (/-S"S,,'S22S"')  . 

Finally.  (A8)  and  (A9)  are  substituted  into  (A4)  and  (A7) 
to  get 

R,A>  =  Sll/M)  +  S'2SU)U2SnSu  ’S21/'  ■" 

+  S'2S{liU2SuliC) 

Ra)  =  S43S</-)(/1S2l/M) 

+  S43S(nl/1S22S(',S34/(°  +  SuIiC). 

S4A,  SiC,  S(A,  and  S(C  can  be  identified  easily  from  the 
above  equations. 
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Microstrip-Fed  Planar  Frequency-Multiplying 

Space  Combiner 

S.  N  \M,  TOMOK.1  UWANO,  member,  ieee,  and  TASUO  ITOH,  fellow,  ilee 


Abstract — A  frequency- multiplying  power  combining  array  has  been 
made  of  slots  in  the  ground  plane  on  a  substrate  with  a  microstrip  fcedline 
on  the  back  side.  The  second  harmonic  generated  by  a  diode  in  each  slot  is 
combined  in  free  space.  A  design  procedure  and  experimental  results  are 
presented. 


1.  INTRODUCTION 

GENERATION  OF  a  coherent  millimeter-wave  signal 
with  a  reasonable  strength  becomes  more  difficult  as 
the  frequency  of  operation  is  increased.  One  of  the  meth¬ 
ods  to  alleviate  this  problem  is  the  use  of  a  power  combin¬ 
ing  technique  [1].  However,  it  is  increasingly  difficult  to 
combine  the  power  outputs  from  more  than  a  few  devices. 
In  many  power  combining  schemes,  very  involved  design 
procedures,  complicated  fabrication,  and  tricky  postfabri¬ 
cation  tuning  are  required.  An  alternative  approach  is  the 
use  of  frequency  multiplication  [2],  In  this  approach,  the 
output  power  level  is  substantially  lower  than  the  input 
signal  level.  In  addition,  the  power  level  that  can  be 
handled  by  a  diode-type  frequency  multiplier  is  quite 
limited.  To  alleviate  many  of  the  deficiencies  described 
above,  a  quasi-optical  frequency-multiplying  space  power 
combiner  has  been  proposed  and  tested  [3],  The  proposed 
structure  contains  in  a  single  component  (1)  a  slot  antenna 
array.  (2)  frequency  multipliers,  and  (3)  a  power  combiner. 
The  structure  has  a  simple  configuration  and  requires  no 
bias  and  no  postproduuion  tuning. 

In  this  paper,  a  new  printed  circuit  configuration  of  the 
frequency-multiplying  space  power  combiner  is  presented. 
In  contrast  to  the  work  in  [3|,  no  waveguide  is  used  for 
feeding.  In  the  new  structure,  the  array  of  slots  on  the 
substrate  is  fed  by  a  miu  ostrip  line  or  a  coplanar  waveguide 
fabricated  on  the  back  side  of  the  substrate.  The  new 
structure  has  other  desirable  features,  such  as  (1)  a  planar 
structure  suitable  for  monolithic  integration,  (2)  large 
coupling,  (3)  controllable  coupling,  and  (4)  flexibility  in 
the  arrangement  of  slots  on  a  planar  circuit. 

Fig.  1  shows  a  typical  configuration  of  the  planar 
multiplier/combiner.  On  the  front  side  of  the  substrate, 
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Fig.  1.  Structure  of  multiplying  slot  array  fed  by  microstrip 

there  is  an  array  of  slots  and  in  each  of  these  a  diode  is 
placed  for  frequency  multiplication.  A  microstrip  line  feed 
on  the  back  side  is  used  in  this  example.  The  incident 
signal  at  the  frequency  f0  is  successively  fed  to  these  slots, 
which  are  typically  a  quarter  wavelength  long.  The  diode 
in  each  slot  generates  a  second  harmonic  of  2/0.  Then  the 
slot  behaves  as  a  half  wave  slot  dipole  and  an  efficient 
radiation  of  a  2/0  signal  takes  place.  The  radiated  powers 
from  all  the  slots  are  combined  in  free  space.  This  paper 
reports  an  experimental  design  of  the  structures  scaled  for 
a  low-frequency  operation  at  the  T-band.  This  scaling 
allows  us  to  use  more  precise  characterizations  than  those 
attainable  at  millimeter-wave  frequencies,  so  that  useful 
design  data  can  be  collected  more  efficiently. 


II.  Design  Procedure 

In  order  to  design  the  proposed  structure,  it  is  necessary 
to  know  the  characteristics  of  the  microstrip-to-slot  tran¬ 
sition.  the  phase  relationship  between  slots,  and  the 
conversion  gain  of  each  slot.  Let  us  first  study  these 
characteristics  in  detail. 


A.  Magnitude  Characteristics  of  \ficrostrip-to-Slut 
Transition 

The  transition  characteristics  have  been  studied  based 
on  the  approximate  equivalent  circuit  shown  in  Fig.  2  (4) 
The  magnitude  characteristics  of  the  coupling  from  the 
microstrip  to  the  slot,  the  insertion  loss  of  the  slot,  and  the 
return  loss  of  the  slot  with  a  50-fi  load  located  at  the 
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n :  1  Z1 

(termination) 

Fig.  2.  Equivalent  circuit  for  one  slot  (/  =  X„8). 


center  of  the  slot  are  given  as  follows: 

Coupling  =  201og|2«Z(cos/Jd  4-  tan  fid  sin  fid ) 

/(( Z0m  +  ZL+Z)(  1  +  jZ0s  tan  fid/Z,))  \ 

0) 

Insertion  loss  =  20 log  |2 Z,  (Z0m  +  Z,  +  Z)\  (2) 

Return  loss  =  20  log  |(  Z0m  -  Z  -  Z, )/( Z0m  +  Z+Z,) 

(3) 

where 

Z  slot  impedance  seen  from  the  microstrip, 

Zu„,  characteristic  impedance  of  the  microstrip, 

Z0,  characteristic  impedance  of  the  slot  line, 

Z,  termination  impedance  (50  ft), 

Z,  slot  impedance  seen  from  the  center  of  the  slot. 

According  to  the  calculated  results,  the  smaller  the  dis¬ 
tance  from  the  center  of  the  slot  to  the  feeding  position, 
the  larger  the  coupling  becomes.  In  order  to  examine  the 
possibility  of  controlling  the  power  coupled  from  the  mi¬ 
crostrip  to  the  slot,  we  measured  the  transition  characteris¬ 
tics.  Measurement  of  microstrip-to-slot  transition  has  been 
carried  out  with  the  conventional  coaxial-to-slotline  transi¬ 
tion  [5]  as  a  probe  for  picking  up  the  signal  over  the  slot  in 
place  of  the  diode.  From  the  measurement  results,  it  is 
found  that  the  control  of  the  power  coupled  from  the 
microstrip  to  the  slot  is  possible  by  controlling  the  distance 
from  the  center  of  the  slot  to  the  feeding  position.  Also,  it 
is  seen  that  the  equivalent  circuit  for  the  transition  gives  a 
prediction  that  agrees  well  with  the  measured  data,  if  n  =  1 
is  chosen,  as  long  as  the  feeding  position  is  not  too  close  to 
the  end  of  the  slot.  The  results  are  shown  in  Fig.  3.  The 
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l'ig  3.  Transition  charactcrisUcs  from  microstrip  to  slot,  (a)  Coupling 
(O.  insertion  loss  ( ),  and  return  loss  ( /., )  (h\  Phase  change  of 
transition  with  load  position. 


dotted  lines  show  the  calculated  results  under  the  assump- 
tion  of  n  =  1  in  Fig.  2. 

B  Phase  Relationship  between  Slots 

The  phase  relationship  between  slots  is  determined  by 
(1 )  the  electrical  length  of  the  microstrip  feed  line  between 
slots,  (2)  the  phase  characteristics  of  the  transition,  and  (3) 
the  square-law  characteristics  of  the  diode.  Among  these 
factors,  the  main  contribution  is  the  electrical  length  of  the 
microstrip  feedline  between  the  slots.  This  quantity  can  be 
controlled  by  the  physical  length  of  the  microstrip  line 
between  the  slots. 

The  phase  characteristics  of  the  transition  can  be 
evaluated  from  the  phase  change  of  the  voltage  from  the 


transition  point  of  the  nucrostrip  to  the  50-ft  load  on  the 
slot  by  the  following  equation: 

Phase  change  =  arg[2«Ztcos/?</  +  tan  fid  sin  fid  ) 

/( (  Z„m  +  Z,  ,  Z  )( 1  +  yZ(),  tan  fid/Z, ))] .  (4) 

The  results  are  plotted  in  Fig.  3(b)  as  a  function  of  the 
position  of  the  load.  From  these  results,  we  find  that  the 
difference  of  the  phase  change  to  the  feeding  position  can 
be  neglected  if  the  load  is  located  at  the  center  of  the  slot. 

A  diode  located  across  the  slot  can  be  thought  to  pick 
up  the  voltage  of  the  input  frequency  signal.  The  voltage 
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on  the  slot  forms  a  standing  wave  pattern  given  by 

V(x)  =  sin  fi{d  +  lh-  x)/sin  p{d  +  lh).  (5) 

Since  this  voltage  has  a  maximum  value  at  the  feeding 
position,  it  is  better  to  locate  the  diode  at  the  feeding 
position  for  n  better  matching  and  a  stronger  coupling.  Ac 
shown  in  Fig.  3(b),  the  phase  change  of  voltage  versus 
feeding  position  of  each  transition  cannot  be  neglected  in 
this  case.  Therefore,  it  is  required  that  the  difference  of  the 
phase  change  at  each  transition  be  considered  for 
determination  of  the  length  of  the  microstrip  line  between 
the  slots. 

In  the  case  of  a  two-element  slot  array,  the  length  of  the 
microstrip  feed  line  between  the  slots  is  Ag/2,  where  \g  is 
the  guide  wavelength  at  the  input  frequency.  However,  the 
diodes  should  be  placed  with  the  same  polarity  at  each  slot 
for  the  broadside  operation.  This  can  be  explained  by  the 
square-law  behavior  of  the  diode  [3], 

C.  Conversion  Gain  of  Each  Slot 

The  conversion  gain  of  a  multiplier  is  defined  as  the 
ratio  of  the  output  power  to  the  input  power.  In  the 
present  case,  the  input  power  is  the  amount  of  power 
coupled  to  a  single  slot  from  the  microstrip  at  the  funda¬ 
mental  frequency.  The  output  power  is  the  amount  of  the 
power  at  the  second  harmonic  frequency,  which  is  radiated 
by  this  slot  into  free  space.  Thus,  we  have  to  evaluate  the 
amounts  of  input  power  and  radiated  power  tc  find  the 
conversion  gain.  The  impedance  of  the  mixer  diode  used  in 
this  paper  is  a  function  of  the  drive  power.  However,  if  the 
diode  is  driven  sufficiently  (more  than  10  dBm),  the 
impedance  is  known  to  be  about  50  fl.  Thus,  the  power 
coupled  into  the  diode  at  the  fundamental  frequency  can 
be  estimated  by  using  C  in  Fig.  3(a).  The  direct  measure¬ 
ment  of  the  amount  of  the  radiated  output  power  is  very 
difficult.  Hence,  we  use  an  indirect  method  based  on  the 
fact  that  the  power  received  at  a  fixed  receiving  antenna  is 
proportional  to  the  output  power  of  a  transmitting  antenna. 
This  relation  is  given  by  [6) 


P,  =  P,  +  G,  +  K 


where 


n 

Pi  ,d 
(S. 4  GHz) 


in. 8  Gil: 


Spectrum 
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Fig.  4  Test  setup  for  measuring  the  multiplication  conversion  gain, 
(a)  DUT  setup,  (b)  Reference  setup. 


Using  (7),  the  power  relation  applied  to  the  DUT  setup 
shown  in  Fig.  4(a)  is  given  by 


Pr.J=  (P,.j  +  cj+Gc)+G,'<j+  K 


where 


Pr  d  receive  power  (dBm)  in  DUT  setup. 

G,  d  antenna  gain  (dB)  of  DUT. 

For  the  reference  antenna,  the  transmit  power  can  be 
expressed  by 

P,.,-P(.,  +  C,  (9) 


where 


P,  receive  power  (  JBm), 

P,  transmit  power  (dBm), 

G,  transmit  antenna  gain  (dB), 

K  constant  determined  by  the  gain  of  receiving  an¬ 
tenna.  distance  between  the  transmitting  antenna 
and  the  receiving  antenna,  and  the  test  frequency. 

For  the  DUT  (device  under  test),  the  transmit  power  can 
be  expressed  by 

P,  ,i  =  P,'j  +  Cj  +  Gt  (7) 

P,  j  transmit  power  (dBm)  from  DUT, 

P,  d  incident  power  (dBm)  to  the  microstrip  feed  line 
in  DUT. 

C(/  power  coupling  factor  (dB)  to  diode  in  DUT, 

G,  conversion  gain  (dB)  of  DUT. 


P,r  transmit  power  (dBm)  of  reference  antenna. 

P'  ,  incident  power  (dBm)  to  the  microstrip  feed  line  in 
reference  antenna, 

Cr  power  coupling  factor  (dB)  to  reference  antenna. 

Using  (9).  the  power  relation  applied  to  reference  setup 
shown  in  Fig.  4(b)  is  given  by 

Prr  =  (P,^  +  Cr)  +  Glr  +  K  (10) 

where 

Pr  r  receive  power  (dBm)  in  reference  setup. 

G,  r  gain  (dB)  of  reference  antenna. 

Using  (8)  and  (10),  the  conversion  gain  can  be  calculated 
from 

G,  =  P,  Pr  r  -  P,  ',  +  P,  ,-Cj  +  C,-G,  /G,  (11) 

D.  Design  Procedure 

Lei  us  consider  the  design  procedure  of  the  multiplying 
slot  array.  Coupling  is  defined  by 

Coupling  =  PJink  -  P ,npul  [dB],  (12) 


v.v.v,  s  v 


1274 


li:i:i.  IRANXAITIONSON  Ml(  KOWAVI:  NIIOKY  ANIi  .i  (  IINIOUI-S.  VOI ..  MIT-35.  NO.  12.  DK  IMHI  K  |VX7 


For  simplicity,  we  locate  all  the  diodes  at  the  center  of 
each  slot  and  design  an  equally  excited,  broadside 
multiplying  slot  array  by  applying  (12)  to  the  measured 
results  in  Fig.  3.  First,  from  the  power  requirement  of  the 
diouc  in  the  first  slot,  we  calculate  the  necessary  coupling 
by  (12);  then  we  evaluate  the  location  of  the  microstrip- 
to-slot  transition  from  C  in  Fig.  3(a).  We  can  then  find  out 
the  amount  of  power  left  over  on  the  microstrip  line  from 
L,  in  Fig.  3(a).  Next,  we  find  the  location  of  the  transition 
of  the  second  slot  so  that  the  power  provided  to  the  second 
diode  is  equal  to  that  for  the  first  diode.  Then,  we  deter¬ 
mine  the  length  of  the  microstrip  to  the  second  slot  to 
obtain  an  appropriate  phase  shift  between  slots.  This 
process  is  repeated  until  all  the  slots  are  taken  care  of. 
Notice  from  Fig.  1  that  the  location  of  the  transition 
moves  inward  as  the  feed  signal  progresses  along  the 
microstrip  line.  The  decreasing  power  available  from  the 
feed  line  is  compensated  for  by  an  increasing  degree  of 
coupling  so  that  all  the  diodes  are  given  an  equal  excitation 
power.  It  is  possible  to  design  various  multiplying  slot 
array  antennas  with  a  suitable  arrangement  of  multiplying 
slots  in  the  ground  plane  with  the  proper  excitation  powers 
and  phase  differences  [6].  These  design  quantities  can  be 
obtained  by  controlling  the  locations  of  the  transitions  and 
the  lengths  of  the  microstrip  between  the  slots. 


III.  Experimental  Results 

For  frequency  multiplication,  nonlinear  reactive  elements 
such  as  varactor  diodes  are  desirable  for  achieving  good 
performance.  However,  only  the  Schottky-barrier  diodes 
designed  for  A'-band  mixers  have  been  available  to  the 
authors.  The  conversion  gain,  Gc,  measured  by  the  method 
described  above  for  the  structure  with  the  ND5051(NEC) 
diode  has  been  -6.2  dB.  In  this  experiment,  the  power 
coupling  factor  to  the  reference  antenna,  Cr,  can  be 
estimated  by 

Cr  =  10 log  [{Pin  -  Pnf  -  P lfa)//> ;„]  (13) 


P,„  incident  power  (mW)  to  microstrip  feed  line  in 
reference  antenna, 

Pwt  reflected  power  (mW)  at  input  port  in  reference 
antenna, 

Py,u  transferred  power  (mW)  to  the  matching  load  in 
reference  antenna. 

An  estimated  power  coupling  factor  of  -7  dB  results. 
Note  that  the  only  difference  in  antenna  structure  between 
the  DUT  and  the  reference  antenna  is  that  there  is  a  diode 
in  the  DUT  while  there  is  no  diode  in  the  reference 
antenna.  This  fact  causes  a  difference  in  the  gain  of  the 
two  antennas.  The  difference  of  the  two  antenna  gains  can 
be  estimated  by  measuring  the  reduction  of  the  receiving 
power  in  the  reference  setup  shown  in  Fig.  4(b)  when  a 
diode  is  located  at  the  center  of  the  slot  of  the  reference 
antenna.  About  10  dB  reduction  in  the  antenna  gain  due  to 
the  loading  effect  of  the  diode  is  measured. 


For  the  other  higher  order  harmonics  generated  by  the 
diode,  the  slot  can  also  be  thought  of  as  a  siot  dipole 
antenna.  However  it  is  observed  that  the  receiving  power 
of  the  third  harmonic  rignal  is  lowei  nan  that  of  the 
second  harmonic  signal  by  more  than  10  dB. 

In  order  to  design  the  multiplying  slot  array  antenna,  it 
is  necessary  to  know  the  characteristics  of  the  element 
radiation  pattern  of  the  single  multiplying  slot  antenna. 
First,  the  effect  of  the  ground  plane  size  is  evaluated.  Fig. 
5(a)  shows  the  radiation  patterns  of  a  single  slot  antenna 
with  a  ground  plane  of  different  sizes.  It  is  found  that  the 
length  of  one  side  of  the  ground  plane  should  be  larger 
than  seven  wavelengths  at  the  radiating  frequency  (about 
20  cm  at  10.8  GHz).  Otherwise,  the  element  radiation 
pattern  of  the  slot  antenna  has  several  peaks  and  valleys. 
This  effect  of  the  finite  ground  plane  on  the  radiation  by  a 
slot  can  be  explained  theoretically  by  the  uniform  theory 
of  diffraction  [6],  Next,  the  effect  of  the  location  of  the 
transition  is  tested.  Fig.  5(b)  shows  the  radiation  patterns 
of  a  single  multiplying  slot  antenna  with  the  transition  at 
different  locations.  We  notice  that  the  position  of  the 
feeding  point  has  a  negligible  influence  on  the  radiation 
pattern.  This  is  due  to  the  fact  that  for  the  second  harmonic 
signal,  the  diode  located  at  the  center  of  the  slot  can  be 
considered  a  source.  Even  though  its  output  power  is 
dependent  on  the  feeding  position,  its  position  on  the  slot 
is  not  changed  with  the  feeding  position.  Therefore,  the 
radiation  power  of  the  multiplying  slot  antenna  is  changed 
with  the  feeding  position  but  the  radiation  pattern  for  the 
front  space  of  the  slot  seems  to  be  insensitive  to  the 
position  of  the  narrow  inicrostnp  feed  line  on  the  back 
side. 

For  the  loading  effect  of  the  diode  on  the  radiation 
pattern,  the  radiation  patterns  of  the  single  slot  antenna 
with  the  diode  and  without  the  diode  are  measured.  By 
comparing  Fig.  5(b)  with  Fig.  5(a),  we  notice  that  rhe 
diode  located  at  the  center  of  the  slot  reduces  the  number 
of  undulations  in  the  radiation  pattern.  However,  the 
overall  radiation  pattern  .s  not  changed. 

Based  on  the  experimental  transition  characteristics  and 
the  elementary  radiation  pattern  of  the  single  slot,  we 
designed  and  tested  several  uniformly  excited  broadside 
multiplying  slot  arrays.  Figs.  5(c)  and  (d)  show  the  radiation 
patterns  of  a  two-element  //-plane  slot  array  antenna  for 
different  spacings  between  the  slots.  It  is  observed  that  the 
radiation  pattern  of  an  //-plane  array  with  more  closely 
spaced  elements  shifts  more  from  the  broadside  direction. 
This  is  due  to  the  increased  internal  mutual  interference 
between  the  slots.  The  spacing  between  //-plane  slots 
should  be  more  than  one  wavelength  at  the  second 
harmonic  frequency  to  avoid  interference.  However,  this 
requirement  produces  grating  lobes  in  pairs  for  a  broad¬ 
side  antenna.  For  the  Jesign  of  an  array  with  closer 
spacing,  it  is  necessary  to  measure  the  element  radiation 
pattern  including  the  effect  of  internal  interference.  Figs. 
5(e)  and  (f)  show  the  radiation  patterns  of  a  two-element 
£-plane  slot  array  antenna  with  different  degrees  of  cou¬ 
pling.  It  is  observed  that  the  radiation  pattern  of  an  array 
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Fig  5.  Measured  radiation  patterns  of  various  multiplying  slot  anten¬ 
nas.  (The  magnitude  is  plotteo  in  linear  scale.)  (a)  Radiation  pattern  of 
a  single  multiplving  slot  antenna  for  different  ground  plane  sizes 
(without  diode),  (b)  Radiato  n  pattern  of  a  single  multiplying  slot 
antenna  for  different  feeding  positions  (with  diode.  20x20).  (e)  Radia¬ 
tion  pattern  of  a  two-element  //-plane  slot  array  antenna  (spacing  - 
A /2)  (d)  Radiation  pattern  of  a  two-element  //-plane  slot  array 
antenna  (spacing  =  A)  (c)  Radiation  pattern  of  a  two-element  E-planc 
slot  array  antenna  for  large  coupling  (  -  5  d B)  (f)  Radiation  pattern  of 
a  two-element  /(-plane  slot  array  antenna  for  small  coupling  (  -  10  dB) 
(g)  Radiation  pattern  of  a  2x4  //-plane  slot  array  antenna  (h) 
Radiation  pattern  of  a  2x4  '  -plane  slot  array  antenna.  (The  dotted 
lines  indicate  calculated  radin"on  pattern.) 


with  more  strongly  excited  slot  elements  shifts  more  from 
the  broadside  direction.  I  his  is  because  the  reflected  power 
at  the  transition  is  increased  as  the  coupling  increases  and 
because  the  neglected  phase  shift  at  the  microstrip-to-slol 
transition  is  increased  the  difference  of  the  feeding 
positions  at  the  two  slots  increases. 

Figs.  5(g)  and  (h)  show  the  radiation  patterns  of  the 
2x4  //-plane  and  the  2  <4  £-plane  slot  arrays  with  the 
calculated  radiation  patterns,  respectively.  The  measured 
radiation  patterns  agree  well  with  the  theoretical  results 
except  for  the  //-plane  radiation  pattern  of  the  2x4 
//-plane  slot  array.  The  retson  for  this  is  explained  above. 
These  arrays  are  excited  with  the  incident  power  of  26 
dBm  at  5.4  GHz.  The  power  fed  to  each  diode  is  estimated 
to  be  11  dBm.  Fig.  6  shows  a  photograph  of  the  2x4 
//-plane  multiplying  slot  array  antenna  for  which  the 
radiation  pattern  is  already  plotted  in  Fig.  5(g). 


(b) 

Fig  6.  Pictures  of  a  2x4  slot  array  antenna,  (ai  Slot  arras  pattern 
(M  Mii'rnvt rin  fYvrl  n.nttcrn 


IV.  Conclusions 

By  using  experimental  data  which  characterize  the  tran¬ 
sition  from  the  microstrip  to  the  slot,  we  designed  a  planar 
nuiltiplier/combiner  in  the  form  of  a  slot  array  fed  by  a 
microstrip  on  the  back  side  of  the  substrate.  The  results 
show  the  feasibility  of  the  proposed  structure.  The  new 
structure  has  desirable  features,  such  as  controllable  power 
coupling  and  flexibility  in  array  geometry  on  a  planar 
substrate,  as  well  as  the  features  of  a  general  multiplying 
slot  array  [3].  The  present  structure  is  suited  for  microwave 
and  millimeter-wave  integrated  circuit  applications. 
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ABSTRACT 

This  paper  proposes  a  new  crosstie  overlay  slow-wave  structure  for 
a  traveling-wave  electro-optical  modulator  and  presents  a  design  method 
including  design  considerations  for  the  optimum  structure.  This  structure 
can  be  designed  to  satisfy  the  phase  velocity  and  impedance  matching 
conditions  simultaneously  over  a  broad  spectral  bandwidth.  This  idea  can 
be  applied  to  a  number  of  electrode  structures  for  broadband 
traveling-wave  optical  modulators. 


INTRODUCTION 

Very  high-speed  optical  modulators  operating  up  to  millimeter  wave 
frequencies  are  required  for  high-bit  rate  optical  communication  and 
high-speed  optical  signal  processing.  The  direct  modulation  of 
semiconductor  lasers  in  that  regime  is  not  effective  due  to  the  cavity 
structure  and  needs  a  short  cavity  to  maintain  the  single  longitudinal  mode. 

*This  work  was  supported  by  the  U.S.  Army  Research  Office  under 
contract  DAAG-29-84-K-0076. 
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External  laser  modulation  using  optical  modulators,  however,  can  be  used 
effectively  at  such  high  frequencies  while  maintaining  the  cavity  mode  and 
high  purity  spectral  characteristics. 

In  particular,  traveling- wave  optical  modulators  are  very  useful  for 
such  external  modulation.  This  is  because  their  modulation  depths  are 
linearly  proportional  to  the  interaction  length  of  the  crystal  when  there  is  no 
phase  mismatch  between  the  modulation  and  light  waves  along  the  crystal. 
In  order  to  maximize  their  modulation  bandwidths,  the  phase  and 
impedance  matching  conditions  must  be  kept  over  a  wideband  and  along 
the  crystal  through  a  proper  design  of  traveling- wave  structures!  1  ]. 

Recently,  a  two-layer  traveling-wave  type  modulator  structure  using 
thin  LiNb03  crystal  on  a  dielectric  subtrate  was  proposed  to  satisfy  the 
above  two  matching  conditions  simultaneously  [2].  Since  the  phase  velocity 
of  the  modulation  wave  is  much  slower  than  that  of  the  light  wave  in  the 
case  of  the  LiNb03  crystal,  the  velocity  of  the  modulation  wave  was 
increased  using  a  substrate  layer  which  has  a  low  dielectric  constant  and 
supports  the  thin  LiNb03  layer.  However,  such  a  two-layer  structure  is 
intrinsically  dispersive  because  the  modulation  fields  are  gradually 
concentrated  into  the  upper  LiNb03  layer  having  a  high  dielectric  constant 
with  the  increase  of  frequency. 

Although  waveguide  modulators  on  semiconductor  substrates  offer 
the  potential  for  monolithic  integration  with  sources  and  electrical  driver 
circuitry,  no  semiconductor  waveguide  modulators  with  traveling-wave 
electrodes  have  been  reported  to  date[3].  This  is  because  the  electro-optic 
coefficients  of  semiconductor  materials  are  relatively  small  compared  to 
LiNb 03.  Their  Schottky  junction  capacitances,  depending  on  the  applied 
voltage,  complicate  the  design  of  broadband  traveling- wave  structures, 
Nevertheless,  the  low  dielectric  constants  and  the  small  velocity  mismatch 
for  both  GaAs  and  InP  are  very  attractive  because  they  allow  a  narrow 
electrode  spacing  for  a  given  charateristic  impedance  and  make  the  phase 
velocity  matching  easy.  However,  conventional  coplanar  traveling-wave 


electrodes  on  semiconductor  materials  cannot  be  designed  to  achieve  this 
phase  matching  because  the  phase  velocity  of  the  modulation  wave  is  in  all 
cases  higher  than  that  of  the  light  wave  due  to  the  field  in  the  air  side. 

This  paper  proposes  an  almost  non-dispersive  slow  wave  structure 
for  a  wideband  traveling-wave  type  optical  modulator.  This  structure  can 
be  designed  to  satisfy  the  phase  and  the  impedance  matching  conditions 
simultaneously. 


CROSSTIE  OVERLAY  SLOW-WAVE  STRUCTURE 
The  proposed  slow-wave  structure  shown  in  Fig.  1  consists  of  an 
infinite  array  of  a  more  capacitive  Section  A  and  a  more  inductive  Section 
B  on  the  coplanar  waveguide[4].  Since  the  period  chosen  is  much  smaller 
than  the  propagation  wavelength,  this  structure  behaves  as  a  uniform 
transmission  line.  In  this  structure,  we  can  increase  the  effective 
capacitance(  C )  and  inductance(  L)  of  this  transmission  line  independently 
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Fig. - 1  Proposed  Slow-Wave  Optical  Modulator  Structure 


by  a  spatial  separation  of  electric  and  magnetic  energies  into  two  different 
sections.  The  shunt  capacitance  and  the  series  inductance  in  the  equivalent 
circuit  of  the  transmission  line  corresponds  to  the  effective  capacitance  and 
inductance  dominating  in  Section  A  and  Section  B,  respectively. 
Therefore,  the  propagation  constant  P ),  which  is  proportional  to  VLC, 
can  be  increased  above  that  of  the  common  coplanar  waveguide.  The 
slow-wave  factor  ( p/p0,where  P0  is  the  free  space  wavelength )  in  the 
range  of  from  2.5  to  12  can  be  obtained  with  an  appropriate  design  of  the 
structure. 

The  numerical  analysis  of  this  structure  has  been  carried  out  by  a 
full-wave  analysis  for  each  section  using  the  standard  spectral  domain 
method[5]  and  by  an  application  of  Floquet's  theorem  to  the  periodic  array. 
The  dielectric  constant  is  assumed  to  be  uniform  in  the  substrate  region  at 
modulating  frequencies  up  to  the  millimeter  wave  frequency  range,  because 
a  small  variation  of  the  dielectric  constant  is  enough  to  form  a  well  guiding 
optical  waveguide.  Although  the  dielectric  loss  is  negligible  for  the 
semi-insulating  GaAs  substrate,  the  conductor  loss  cannot  be  neglected  in 
this  small  device  structure.  The  real  and  imaginary  parts  of  the  complex 
propagation  constants  corresponding  to  the  phase  and  the  attenuation 
constants  can  be  calculated  using  the  spectral  domain  analysis  and  the 
incremental  inductance  formula[6j,  respectively. 

While  the  calculated  complex  propagation  constants(  PA,PB )  of 
Section  A  and  Section  B  are  close  to  each  other  and  below  that  of  the  light 
wave  in  the  channel  waveguide,  the  overall  propagation  constant  p )  of  the 
periodic  array  that  was  calculated  using  the  Floquet’s  theorem  can  be 
increased  by  the  different  characteristic  impedances  of  the  two  sections. 
The  following  dispersion  and  impedance  equations  are  derived  by  applying 
the  Floquet's  theorem  and  continuity  conditions  of  voltage  and  current  at 
each  period  junction.  The  discontinuity  effects  at  the  junction  are  assumed 
to  be  small  for  the  simplification. 


cosCyl)  =  cos(yAIA)  cos(yr1b)  -  (1/2)(ZA/ZB+ZR/ZA)  sin(YAlA)  sin(yBlB) 
z  =  F  ZA  sin(YAlA)  cos(Yb!b)  +  Zb  C0S(Ya1a)  sin(YBlB)  1  1/2 
l^l/ZA)  sin(YAlA)  cos(yb1b)  +  (1/ZB)  cos(ya1a)  sin(YB!;!)^ 

where  %  Z  =  overall  complex  propagation  constant  and  characteristic 
impedance 

Ya,  Yb  -  complex  propagation  constant  of  Sections  A  and  B, 
respectively 

ZA.  ZB  =  characteristic  impedances  of  Sections  A  and  B,  respectively 

1  =  1A+1B,  with  1A  and  1B  being  the  lengths  of  Sections  A  and  B, 
respectively 

Assuming  for  simplicity  that  yAlA,  YB1B  «  1,  YA  =  Yb  and  ZA  «  Zp  the 
above  equations  can  be  reduced  to  the  following  simple  expressions 

(  1+T)p  -  pA(  1+TK  +  T2)  1/2 

(  1+T)2a  -  (PA/P)[aA  +  T2aB  +  KT(aA  +  aB)/2] 

Z  -  [ZAZgTK/(K+T))1/2 

Since  p  and  Z  are  approximately  proportional  to  VzB/ZA  and  VZAZB 
respectively,  in  the  event  of  1  «  ZB/ZA  and  1A=  1B,  we  can  easily  increase 
the  propagation  constant  while  maintaining  a  given  characteristic 
impedance.  Therefore,  this  property  can  be  applied  very  effectively  in 
order  to  achieve  the  simultaneous  matching  of  the  pha.^e  and  the 
impedance,  that  is  required  in  the  optimum  electrode  design  of 
traveling-wave  optical  modulators.  Furthermore,  since  this  slow-wave 
structure  is  almost  dispersionless  due  to  the  uniform  transmission  line 
behavior,  the  structure  is  very  useful  for  the  design  of  wideband 
traveling-wave  optical  modulators.  This  slow-wave  phenomenon  can  also 
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be  applied  to  a  number  of  traveling-wave  electrode  structures  for  the 
optimum  design  of  a  wideband  optical  modulator. 


PHASE  MATCHED  OPTICAL  MODULATOR 
The  proposed  slow-wave  optical  modulator  can  be  fabricated  as 
follows  as  shown  in  Fig.l.  An  AlGaAs  layer  is  grown  on  a 
semi-insulating  GaAs  substrate.  GaAs  optical  channel  waveguides  are 
formed  on  the  AlGaAs  layer  by  a  buried  structure  technique[7].  Coplanar 
waveguide  electrodes  are  fabricated  on  the  top  surface  of  the  AlGaAs  layer. 
Next,  an  insulating  dielectric  layer  is  deposited.  Finally,  periodic  metal 
crosstie  overlays  are  fabricated  on  the  dielectric  layer.  Therefore,  the 
dielectric  loss  of  the  modulation  wave  and  the  optical  loss  due  to  free 
carrier  absorption  can  be  avoided  in  this  structure.  The  localized  effective 
index  change(8n~n3rE/2),  which  is  formed  by  the  electro-optic  effect 
for  the  traveling  electric  field(  E  )  of  the  modulation  wave,  generates  the 
local  phase  shift  of  the  optical  wave  over  the  interaction  length(  L )[  1].  In 
order  to  accumulate  the  local  phase  shift  along  the  waveguide,  the  phase 
matching  between  the  optical  and  modulation  waves  can  be  accomplished 
by  our  proposed  slow-wave  structure. 


Design  of  Slow-Wave  Coplanar  Electrode 
This  slow-wave  structure  has  many  design  parameters  including  a 
coplanar  waveguide,  overlay  crossties  and  a  dielectric  layer.  This  allowes  a 
wide  design  flexibility  for  the  simultaneous  matching  of  the  phase  velocity 
and  the  impedance. 

The  two  design  requirements  for  the  simultaneous  matching  of  the 
phase  velocity  and  the  impedance  are  given  by 

Z  =  50  (Cl) 
nmc  =  n'e 


where  nme  and  n^are  the  effective  index  of  the  modulation  and  light 
waves,  respectively. 

Additional  design  considerations^]  should  be  included  for  high 
modulation  efficiency,  low  conductor  loss  and  realizable  dimensions.  High 
field  intensity  can  be  achieved  by  the  minimum  realizable  electrode  gap 
while  low  current  density  requires  thick  and  wide  electrodes.  Since  the 
applied  electric  fields  in  sections  A  and  B  are  mostly  concentrated  in 
dielectric  and  substrate  layers,  respectively,  it  is  necessary  to  increase  the 
length  ratio  between  the  two  sections(lB/lA)  in  order  to  achieve  the  high 
effective  overlap  between  the  applied  field  and  the  optical  mode.  The  one 
period  length(lA+lB)  for  a  fixed  ratio(  1B/1A)  affects  only  the  bandwidth 
and  is  limited  by  the  minimum  achievable  feature  size  of  electrodes.  Silicon 
Nitride  is  chosen  for  a  dielectric  material  because  of  its  high  dielectric 
constant  and  high  field  strength  which  allow  high  slow-wave  effect  and 
high  field  operation,  respectively. 


Design  of  Optical  Channel  Waveguide 
The  optical  waveguide  must  be  designed  interactively  with  the 
electrode  design  to  enhance  the  field  efficiency  of  the  modulation  wave. 
The  main  factors  for  the  optimum  design  are  the  effective  index,  the  mode 
size  and  the  relatve  alignment  to  the  electrodes. 

Since  the  composition  ratio(x)  of  the  Al^Ga,  xAs  layer  can  be 
accurately  controlled  over  a  wide  range  using  epitaxial  growing  techniques, 
the  refractive  index  difference(  An  «  0.62  Ax,  where  Ax  is  the  composition 
ratio  difference)  between  GaAs  channel  and  AlGaAs  clad  of  the 
waveguide  can  be  made  large  enough  for  the  high  field  confinement  of  the 
optical  mode[8].  The  mode  penetration  into  the  dielectric  layer  is  also  quite 
small  due  to  the  high  refractive  index  difference  between  the  dielectric  and 
the  GaAs  layers.  Therefore,  the  effective  index(  n'c)  and  the  horizontal 
and  the  vertical  mode  sizes(  wh,wv )  can  be  simply  expressed  as  shown 


below  for  the  well  guided  fundamental  mode  in  rectangular  waveguides[9]. 
n'c  =  [n,2-(X/2a  )2-(X/2b  )2]1/2 

wh  =  a  +  0.23  X  /  ( An-n2 
wv  =  b  +  0. 12  X  /  (  An-n2  )1/^ 

where  n^,  n2  =  refractive  index  of  GaAs  channel  and  AIGaAs  clad 
.respectively 

a,  b  =  horizontal  and  vertical  length  of  the  rectangular  waveguide 
.respectively 

X  =  free  space  wavelength  of  the  optical  wave 

Although  small  mode  size  and  low  effective  index  are  advantageous 
for  high  field  overlap  and  easy  phase  matching,  the  minimum  mode  size 
should  be  comparable  with  those  of  typical  single  mode  fibers  ( »  7  ^im ) 
for  low  fiber  coupling  loss.  Therefore,  in  the  optimum  design  there  is  a 
tradeoff  between  the  high  field  overlap  and  the  low  coupling  loss.  Since 
the  minimum  electrode  gap  in  this  structure  is  limited  by  the  impedance 
matching  to  10  tim,  the  practical  ratio  between  the  mode  size  and  the 
electrode  gap  lies  around  0.7.  Therefore,  the  waveguide  placed  between 
the  electrodes,  which  uses  the  horizontal  component  of  the  modulation 
field,  is  effective  for  high  field  overlap(  T  =0.6)  and  has  low  optical 
loss[10]. 


NUMERICAL  RESULTS  AND  DISCUSSION 
There  are  many  possible  combinations  of  multiple  design  parameters 
for  the  simultaneous  matching  of  the  phase  velocity  and  impedance.  One 
combination  listed  in  Table- 1  is  obtained  through  the  above  mentioned 
design  considerations  and  the  numerical  calculations. 


Table-1  Design  Parameters  lor  the  Slow-Wave  Modulator 
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Table-2  Numerical  Results  for  the  Design  of  Table-1 

(at  1  GHz  modulation  frequency  and  1.15  )im  light  wavelength) 
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The  effective  index  dispersion  and  the  characteristic  impedance 
change  for  the  design  in  Table- 1  are  shown  in  Fig-2.  From  this  figure  it 
can  be  seen  that  this  slow-wave  optical  modulator  has  a  very  wide 
bandwidth  for  the  simultaneous  phase  and  impedance  matching.  In  order  to 
see  the  frequency  characteristics  more  clearly,  we  refer  to  the 
phase-modulation  reduction  factorfll]  given  by 

r\  =  sin  ( u )  /  u 
u  =  :t  f  L  (  nmc(f)  -  n'c )  /  c 

wh^re  f,  L  and  c  are  the  modulation  frequency,  the  modulator  length  and 

the  free  space  velocity,  respectively.  In  Fig-3  the  phase-modulation 
reduction  factor  ( rj )  is  plotted  against  the  frequency-length  product 
f(GHz)-L(cm).  The  frequency-length  product  at  the  3  dB 
phase-modulation  factor  is  about  320  (GHz-cm),  which  is  about  50  times 
larger  than  those  of  normal  LiNb03  modulators  and  7  times  larger  than 
that  of  a  phase  matched  LiNb03  modulator[2].  Therefore,  this  modulator 
can  be  operated  up  to  the  bandwidth  of  320  GHz  for  the  modulator  length 
of  1  cm. 


Fig.-  2  Dispersion  and  Impedance  Characteristics 
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Conductor  Loss  and  Modulation  Voltage 

The  overall  conductor  loss  per  unit  length  of  this  structure  is  about 
two  times  larger  than  that  of  the  common  CPW  due  to  the  slow-wave 
factor  and  the  high  conductor  loss  in  cross-tie  section.  Therefore,  the 
effective  interaction  length  ( Leff )  is  reduced  from  the  actual  length(  L )  as 
below  in  the  phase  matched  case. 

Lcff  =  (l-e-VO/cXy 
where  av  is  the  voltage  attenuation  constant. 

The  maximum  effective  length  for  the  infinitely  long  modulator(  L  = 
oo)  is  the  inverse  of  the  attenuation  constant.  Therefore,  the  minimum 
voltage  required  for  the  half-wave  modulation(  A(3  =  n  )  of  this  modulator 
can  be  given  by  the  following  equation  in  the  case  of  perfect  phase 
matching. 

v"mir  -  «v/n3rr 

Although  the  theoretical  V  rtmin(  =  10  V  at  1GHz  )  is  a  litttle  larger 
than  those  of  the  LiNb03  modulators  having  assymmetric  coplanar 
electrodes[12],  it  can  be  reduced  by  using  the  assymmetric  coplanar 


slow-wave  structure.  That  allows  a  smaller  electrode  gap  for  a  given 
characteristic  impedance. 


CONCLUSIONS 

In  this  paper,  a  new  crosstie  overlay  slow-wave  structure  is 
proposed  for  wideband  traveling-wave  type  optical  modulators.  The 
simultaneous  matching  of  phase  veloity  and  impedance  has  been 
successfully  achieved  through  the  numerical  analysis  of  the  slow-wave 
electrode  structure  and  the  design  of  the  optical  modulator.  This 
simultaneous  matching  can  be  accomplished  over  a  wide  bandwidth  due  to 
the  uniform  transmission  line  behavior  of  the  proposed  electrode  structure. 


REFERENCES 

[1]  R.  C.  Chen,  Proceedings  of  the  IEEE,  Vol.  58,  pp.  1440-1457, 
October  1970. 

[2]  K.  Atsuki  and  E.  Yamashita,  J.  of  Lightwave  Tech.,  Vol.  LT-5,  pp. 
316-319,  March  1987. 

[3]  R.  Alferness,  IEEE  Trans.  Microwave  Theory  Tech.,  Vol.  MTT-30, 
pp.  1121-1137,  August  1982. 

[4]  T.  H.  Wang  and  T.  Itoh,  1987  IEEE  International  Microwave 
Symposium  Digest,  Vol.  1,  pp.  315-318,  June  1987. 

15]  T.  Itoh,  IEEE  Trans.  Microwave  Theory  Tech.,  Vol.  MTT-28, 
pp.  733-736,  July  1980. 

[6]  R.  A.  Pucel  et  al.,  IEEE  Trans.  Microwave  Theory  Tech.,  Vol. 
MTT-16,  pp.  342-350,  1968. 

[7]  K.  Saito  and  R.  Ito,  IEEE  J.  of  Quantum  Electr.,  QE-16,  pp. 
205-215,  February  1980. 

[8]  S.  Adachi,  J.  Appl.  Phys.,  Vol. 58(3),  pp.  R1-R29,  1  August  1985. 


[9]  E.  A.  J.  Marcatih,  Bell  Sys.  Tech.  /.,  pp.  2071-2103,  September 
1969. 

[10]  D.  Marcuse,  IEEE  J.  of  Quantum  Electr .,  QE-18,  pp. 
393-398,  March  1982. 

[11]  L.  P.  Kaminow  and  J.  Liu,  Proceedings  of  the  IEEE ,  Vol.  51,  pp. 
132-136, January  1963. 

[12]  R.  C.  Alferness  et  al.,  Electron.  Lett.,  Vol.  13,  pp.  309-310,  March 
1986. 


-  '  v  '  ‘ ^  -  “j.  vV  '«*  *  V  Y  V  V  ">%v>  “j"  -  “  >  Y%y  “  .’VV-  V  ' 


VjVj’V- 


